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Abfjcrnct 

Our earlier measurements of the high degree and order (n,m>12) 
gravity in the central Pacific using the method of £acellite--jto-£atei],ite 
^tracking (SST) have been extended with an additional 50 passes of data. 

The SST method utilizes line of sight Doppler tracking between the very 
high geosynchronous ATS-6 spacecraft and the much lower (840 km) orbiting 
GEOS-3 spacecraft. The observed changes in velocity with time are re- 
duced relative to the well-determined low degree and order (n,m^l2) GEM-/0 
field model and accelerations are' found by analytical differentiation of 
the range rates. This new map is essentially identical to the first 
map and we have produced a composite map by combining all 90 passes 
of SST data. The resolution of the map is at worst about 5° and much 
better in most places. A comparison of this map with conventional GEM 
models shows very good agreement. A reduction of the SEASAT altimeter 
data has also been carried out for an additional comparison. Although 
the SEASAT geoid contains much more high frequency information, it agrees 
very well with both the SST and GEM fields. The maps are dominated 
(especially in the east) by a pattern of roughly east-west anomalies 
with a transverse wavelength of about 2000 km. A further comparison 
with regional bathymetric data shows a remarkably close correlation with 
plate age. Each anomaly band is framed by those major fracture zones 
having large offsets and the regular spacing (=10°) of these fractures 
seems to account for the fabric in the gravity field. There are other 
anomalies that are accounted for by hot spots and altogether the immediate 
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fjoucce of nt lc*ast pni*t of chese Pnetfic anomalies Is in Che lithosphere 

icself. It therefore seems chat most of the anomalies in the east half 

of the Pacific could be partly caused simply by regional differences in 

plate age. The amplitude of these geoid or gravity anomalies caused by 

age differences should decrease with absolute plate age, and large anomalies 

('^'3 m) over old, smooth sea floor may indicate a further deeper source 

within or perhaps below the lithosphere. We have a].so considered the possible 

plume size and ascent velocity necessary to supply deep mantle material 

to the upper mantle without complete thermal equilibration. A plume 

4 

emanating from a buoyant layer 100 km thick and 10 times less viscous than 
the surrounding mantle will have a diameter of about 400 km. And it must 
ascend at about 10 cm/yr to arrive still anomalously hot in the uppermost 


mantle. 
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Introduction 


In an earlier prper we useu the method pioneered by W. L. Sjogren 

t 

of sateHlte-to~satellite_,tracking (SST) to measure the gravity field 
over the Pacific (Marsh et al. , 1981). That studyused 40 passes of 
ATS-6/GEOS-3 data to construct a high degree and order (n,m > 111) gravity map. The 

present study includes data from about 50 additional passes to produce 
a more detailed gravity map. 

In brief, this method measures the earth’s gravity field by 
Doppler tracking of a low orbiting satellite (GEOS-3 at 840 km) by a much 
higher (40,000 km) geosynchronous satellite (ATS-6) . Doppler tracking 
furnishes the speed of the low satellite as a function of time or position., These 
range-rates are then converted through differentiation to line-of-sight accelera- 
tions or gravity anomalies. Point values of these anomalies measured along each 

track are then contoured into a map of gravity anomalies « The long wave- 

» 

length components of the SST measurements are removed using a low degree and 
order (n,m£l2) field model, which is very well determined. The residual 
anomalies thus represent harmonics of n,m>12« 

The great advantage of the SST method is that it is simple. No large 
arrays of spherical harmonics need be found or manipulated, and there is 
little chance for the data to become numerically adulterated during pro- 
cessing. The method is limited by the height of the low-orbiting satellite 
and thus is sensitive to wavelengths longer than about 1500 km (ige, n,m jc 24) . 
But as shall be shown, it is this range of the gravity field (12 n,m _< 25) 
that is geophysically particularly interesting. A disadvantage of the method 
is that the map applies only at the altitude of 840 km and must be separately 
downward-continued to give anomalies on the earth's surface ^ This method 
also does not strictly measure the radial component of gravity, but rather 
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a changing ilne-of-sighc component of the field between ATS-6 and <UiOH~3, 

Only within 30-40 degrees of the high satellite's subsateiiite point can these 
an.)(iialie.s be. taken as good measures of the radial field, liven with these limita- 
tions we found in our previous study that with only 15 passes of data the map 
derived of the gravity field over North America is essentially identical to that 
produced by any previous field model, which Involves data from 30 or more 
satellites in addition to surface data. 

The gravity data is not only geophysically useful, but it provides an 
excellent opportunity to check one gravity field against another. For example, in 
our previous study we compared the more conventional Goddard Earth Model (GEMIOB; 
Lerch et al., 1981) and the 8ST gravity fields. We also show for the first time 
comparisons with the SEASAT altimeter-derived geoid over the Pacific, which is 
highly accurate and useful for comparisons and geophysical interpretations. We 
also show that there is often a good correlation between relatively young sea 
floor and gravity and geoid anomalies, the major fracture zones of the Pacific, 
because of their offsets, give the distinctive east-west fabric to the gravity and 
geoid anomalies. Other anomalies correlate with the traces of hot spots, and 
overall it is the interference of these two sources of anomalies that gives rise 
to the distinctive gravity field of the Pacific area. 

SST Map 

A map of satellite tracks used here and in our previous study are shown . 
combined as Figure 1. The coverage is particularly dense over the central 
Pacific. The gravity is determined at 70 km intervals along each track and 
these values have been contoured into a map sho^m as Figure 2. This map shows 
the same assemblage of anomalies as we previously noted, but now there 
is more detail. In particular: The positive Line Islands anomaly, in the 
vicinity of the subsatellite point, is more rectangular and strikes 
nearly east-west. The positive tail on the southeast corner of tliis anomaly 
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is now siv.alLor aiul cssentiaUy Insignificant fri,im •/.(.•ro. The positive 


c 


anomaiy around Hawtiii (see plotted islands) and nortlieast of Hawaii is much 
the same as before except the contours within the anomaly show more detail. 
There is also now a stronger bridge ot positive values between the Cook-Austral 
and Pitcairn anomalies. Otherwise the map is essentially Che same as the previous 
SST map of Marsh ot al,, (1981). 


Comparison with SEASAT Geoid 

Since the geometry of an essentially radial gravity field is the 
same as its geoid, we show for comparison with the SST map a geoid over the same 
area (Figure 3) . This geoid (mean sea surface) was computed using the radar alti- 
meter from the SEASAT satellite and it represents 70 days of data taken between 
July and October 1978. The long wavelength field has again been removed 
by subtracting a low degree and order (n,m ^ 12) field model (PGS-S4; 

Lerch et al. , 1982) from the original data. Each data point used in con- 
structing this map represents an average of 1000 radar measurements taken 
Over a time of 1 sec or a ground distance of 7 km, and each point has a 
precision of 3-4 cm (Tapley et al* , 1982). The accuracy of Che altimeter- 
derived mean sea surface with respect to the center of mass of the earth 
is dominated by radial orbit error. The recent SEASAT ephemerides computed 
at GSFC using laser, and Unified S-Band tracking datr, and the most accurate 
earth gravity and geodetic models have a global rms radial accuracy of 70 cm. 
This set of data therefore contains far more high frequency information over 
the ocean than any previous GEM or SST gravity field. Through gridding and . 
contouring, however, these data have- been smoothed somewU :it by using all 
points within a cap with a radius of five degrees. More details of the 
procedure of computing this mean sea surface are described by Marsh and 


Martin (1982). 
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Although f,he SEASAT geoid contains much more Infonnation than the 
f^5T liuip, tin; tv-u iiiapH sliow eSHentlaUy t(u* same aiUKiaUey, Lc should Ik* 
recalled, however, that the relative zero-level between the two maps may 
be different and some features may appear larger or smaller on one or the 
ocher of the maps. There are nevertheless some significant differencee. 

The Hawaiian anomalies (200°, 20‘’N; 215®, 30“N) are hardly connected on Che 
SEASAT geoid. There is also a fairly strong geoid anomaly all along the chain 
of Hawaiian Islands c The SEASAT geoid shows almost a connection between 
the Hawaiian anomaly and Che Line Islands anomaly to the south, Other 
than a small break near 190® the Line Island geoid anomaly nearly connects 
with the Gilbert-Marshall anomaly and forms a positive band nearly traversing 
Che entire Pacific. On both maps this anomaly forms the backbone of the 
pattern of approximately east-west features with a transverse wavelength of 
about 2000 km. This had been seen in our earlier work but is now particu- 
larly clear in these maps. 

In sum, the SST and SEASAT ',Taps show Che same basic areas of positive 
and negative anomalies and overall correlate closely. The SEASAT geoid 
provides a more complete coverage and over Che Pacific as a whole is 
expected to be more accurate than the SST map. 

Correlations with Fracture Zones and Hot Spot Traces 

We noted in our recent (1981) paper that some of these gravity anomalies 
correlate well with residual depth anomalies. The Hawaiian and Cook-Austral 
anomalies, for example, correlate to various degrees with residual depth 
anomalies delineated by, respectively. Watts (1976) and Crough (1978). 
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And we f?ound through n olmple model that shallower-thnn-tisual a«.eas of 
the sea floor can certainly cause the requisite SST gravity anomaly 
if they are compensated at a depth of about 30-100 km. Some positive 
anomalies (e.g. northeast of Hawaii near 220“ and 30“ ) apparently do not 
correlate with residual depth anomalies. The orderly pattern of east- 
west v*;rending negative anomalies just south of Hawaii and the Line Islands 
anomalies, which nearly span the Pacific, apparently do not correlate with 

negative residual depth anomalies. Extensive residual depth data, however, is 
largely wanting. 

Across fracture zones with large offsets in lithosphere age (=^10 m. y.) 
there is a significant change in thickness of the lithosphere, which is 
reflected by a sudden change in ocean depth. This step change in topography 
is well known to produce a distinctive gravity or geoid anomaly. 

The change with age of these anomalies has been used to study 
the aging and growth of the lithosphere (e.g. Detrick, 1981), and 
these anomalies may persist for great distances. In their study of the 
Mendocino fracture zone, for example, Sandwell and Schubert (1982) showed 
that these effects may persist to ages of 135 m.y., which corresponds to a 
distance near the bend in the Emperor-Hawaiian chain of seamounts. All of 
these studies, however, have used a local array of tracks to study a single 
anomaly or offset. Since the Pacific plate is cut by a regular- pattern 
of fracture zones with large offsets, these may produce a certain fabric in 
the gravity and geoid fields. 

To investigate such correlacions betwee'' the gravity and geoid fields 
and regional bathymetry we first made a map of the 1“ x 1“ depths as 
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.supplied by NOAA (1080) for Che north Pacific. A c*.omp<« risen of Ftj^ure 2 
and 3 with this map shows that major fracture xones essentially frame 
the east-west pattern of geoid and gravity anomalies. The major positive 
geoid anomaly northeast of Hawaii cor'^alates with the relatively young 
lithosphere bounded between the Murray and Molokai fracture zones. Because 
the ocean floor magnetic anomalies closely record the age of the plate 
and thus also its bathymetry, there* should also be a positive correlation 
between geoid anomalies and the age of the ocean basin as deduced from the 
magnetic anomalies of the sea floor. The superposition of geoid on plate 
age is shown as Figure 4. From this combined map there is a clear 
correlation between the gravity anomalies and the age of the lithosphere. 

The large fracture zones form a frame work around the gravity anomalies. 

These correlations suggest that gravity and geoid anomalies not 
associated with active volcanism (i.e. so-called hot spots), especially in 

ft 

the eastern Pacific, may result from differences in age of the sea floor. 

The particular pattern of fracture zones where major offsets occur at 
approximately equally-spaced intervals («1000 km) produces a fabric in the 
geoid that dominates the entire Pacific. 

If it were not for the seamount chains, the gravity and geoid 
fields of the Pacific would appear as approximately east-west bands of 
alternating sign. The topography and associated geoid fields of seamount 
chains are superimposed on this regional plate fabric, which is set at the 
ridges. Although the so-called hot spot seamount chains form at some 
distance from the ridge, other chains (e.g. Line Islands) have apparently 
formed near the ridge (Watts et al., 1980). Thus tho Pacific gravity 
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and geoid fields are principally the result of these two distinctive effects. But 
because the course of the l‘actfic plate for the iii.> . iuis riot been purailel to 
the pattern of fracture ;?one3, the twi,, sets of anomalies interfere. Hence 

we turn to a more detailed examination of the principal gravity and geold 
anomalies and their relation to fracture zones, plate age, and hot spot activity. 

Hawaiian Anomaly ; The Hawaiian chain of seamounts, for example, cuts 
across the region bounded by the Murray and Molokai fracture zones which other- 
wise would be a positive band of gravity, and also across the negative region 
to the north between the Mendocino and Murray fracture zones (see Figure 5). 

This may have the resultant effect of producing the relative low in the geoid 
along the Hawaiian seamounts near the Midway Islands. This could have important 
implication* for understanding the correlation of residual depth anomalies and 
geoid anomalies. 

We found in our earlier (1981) study that the present gravity and geoid 

ft 

fields shox/ positive anomalies at the bend in the Hawaiian-Emperor seamount 
chain. We also suggested that this anomaly may be an obstacle to Crough’s 
(1978) model of lithospheric thinning near Hawaii, through local reheating, 
followed by a monotonic decay with cooling of both tope ^raphy and geoid 
anomaly. If this anomaly at the bend is actually a continuation of 
the main Hawaiian chain anomaly, then the orginal Hawaiian anomaly may take 
longer than expected to decay. The eastern portion of this anomaly 
correlates with the shallow sea floor of the Hess Rise (e.g. Nemoto and 
Kroenke, 1981) but the western portion does not. Following Crough's model , 
however, even an original irregularity in the sea floor should be erased 
by the reheating and buoyant rise of the lithosphere to a virtual age of about 
25 m.y. It is therefore not clear whether it is the positive anomaly at 
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the bend or the necking of the lluwation Anomaly near Midway that ts 
anomalous. For this model this differentiation could be important. 

Line Islands Anomaly ; A' similar interference occurs between the 
east-west trending, positive Line Islands geold anomaly band (Marsh 
et al., 1981) south of Hawaii and the trend of the Line Islands them- 
selves. Watts et al. (1980) show that the gravity over the Line Islands 
is consistent with their geologic age of about 97 m.y., implying that 
they formed at a ridge. The continuation to the southeast of the trend 
of the Line Islands, however, has been volcanically active in the last 
1.3 m.y. (Duncan and McDougcll, 1974). Residual depth anomalies have been 
computed for this area by Crough and Jarrard (1981), who found that this anomaly 
decays roughly as predicted for a hot-spot trace that includes the incon- 
sistently older Line Islands, This may imply that the Line Islands have 
been reheated, perhaps as a result of the more recent volcanism of the 
Cross Trend Line or from much earlier, unnoticed , intrusive activity of 
the Marquesas Islands. 

This residual depth anomaly, the Marquesas-Line Swell (Crough and 
Jarrard, 1981), strikes northwesterly but more northerly than that of the 
trace of Pacific plate motion for the last 40 m.y. The associated geoid 
anomaly shown by Crough and Jarrard was derived from the GEOS-3 geoid 
of Brace (1977) by subtracting a regional geoid described by a seventh 
degree polynomial. This geoid anomaly also shows a IWW trend, a trend 
which is not evident in the more precise SEASAT geoid of Figure 3 nor in the 
SST and GEM gravity fields, each of -which have had a twelfth degree global 
field removed. The trend of these latter fields is markedly east-west 
which persists for about AOOO km. In fact, a principal A m anomaly of 
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with tha Mnrqueyas depth anomaly, is absent in che SKASAT data. This 
descrcpancy was described in our earlier (1981) 8ST study and was attributed 
to the sparseness of GEOS-3 altimeter data in this region . The extensive 
SEASAT data confirms this interpretation. 

The WNW trend associated with the Line Islands is evident on the 
northwest corner of this regional anomaly, but any correlation with the 
Marquesas-Line Swell itself is less clear. The choice of removal of the 
regional field is evidently crucial in producing a geoid anomaly that 
correlates with this residual depth anomaly. 

Regardleiss of the age of the Line Islands, they do produce a geoid 

anomaly trending NW, which might be expected, on grounds of seamount 

frequency, to continue through to the Hawaiian anomaly. These anomalies 

nearly do connect, but they are separated by a distinctive regional 

* 

negative anomaly that spans the entire Pacific. This negative band has 
its origin on the east between the Clipperton and Murray fracture zones, 
which together frame relatively older lithosphere. The anomaly does not 
carefully follow these bounds, which is probably due to interference with 
the geoid anomaly of the Line Islands, and it cuts across moderately shallow 
(<5000 m) sea floor' near 20“ and 170-180“. 

Co ok- Austral Anomaly : This gravity and geoid anomaly centered near 

-20“,205“ correlates roughly with the residual depth anomaly computed here 
by Crough (1978). On close inspection, the principal part of the depth 
anomaly (-800 m) is near -28“,215“ and trends northwesterly , which correlates 
with a tail of the geoid anomaly to the southeast. The amplitude of this 
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significantly less than might be expected judging from past coccelatlonu. 

The much more regional nature of the geoid ano gravity anomalies may 
in part be due to the proximity of the Tonga-Kernutdec trench and its 
associated outer rise. It is also clear in this region that there Is 
no positive anomaly over the Manihiki Plateau, which is just northeast 
of Samoa. 

Early radiometric dating of lavas of the Cook-Austral Islands 
showed a broad consistency with the hot spot hypothesis, but four of these 
islands were updated and three others had only dates of unlocated 
samples. Turner and Jarrard (1982) have recently reported ninety four 
additional dates that supply this missing information. The ages of this 
t,'fii»it/ now appear to be altogether inconsistent with those predicted 
by a simple hot spot migration. They instead suggest that the results 

ft 

give stronger support to a "hot line" hypotheses, which may stretch from 
the Samoan chain through the Cook-Australs and on to Pitcairn Island. A 
"hot line" could relate these volcanisms, but there seems to be little 
evidence of it in the gravity or geoid of this area. 

Pitcairn Anomaly : This anomaly is bounded on the east by the east 

Pacific rise and oft the north and south by two fracture zones with large 
offsets. It does not form a continuous band westward to the Cook-Austral 
anomaly, but is bounded on the west by a shallow negative anomaly (see more below) . 

Marshall-Gilbert Anomaly : This diffuse positive anomaly at 0",180“ 

is separated by a distinct but shallow low from the Line Islands anomaly 
to the east. The sea floor in this region is not highly populated with 
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seamouats and is generally 5000-6000 m deep. The Cecconic history of 
this region oL thr i'aciCic l,« nut clear, but rhe dlsitnctly linear, 
north-northwest trending facing edges of these anomalies (i.e. Gilbert- 
Marshall and Line Islands) may suggest the presence of NNIV trending 
fracture zones or other plate fabric. Some confirmation of this comes from 
the correlations of magnetic anomalies in this region whose strike is 
normal to the strike of the inferred fracture zones (e.g. Pitman et al. , 

1974). The chain of islands made up of the Marshall and Gilbert Islands, 
although also striking in this (NNW) direction, does not have an associated 
distinctive geoid anomaly. This general positive anomaly may therefore 
be caused by a region of relatively young sea floor (see below) . 

Geoid and Gravity Anomalies Between Fracture Zones 

The step-like change in geoid height across fracture zones separating 
plates of different ages has been ’used to study the aging of the lithosphere 
by Crough (1979), Detrick (1981), and Sandwell and Schubert (1982). The 
amplitude of these anomalies is proportional to the offset in age, and the 
history of the anomaly can be used to evaluate the cooling characteristics 
and growth of the lithosphere. These results suggest deviations beginning 
at ages as early as 20-40 m.y. from the model of the lithosphere as a thermal 
boundary layer (Sandwell and Schubert, 1982). l>flien modeled instead by the 
cooling of a plate of constant thickness, as suggested by Parsons and 
Sclater (1977) , the closest fit to the change in slope of the geoid anomaly 
with age is for a plate thickness of 125 km (Sandwell and Schubert, 19&2). 

Although in these studies the proper separation of the anomaly due to 
the fracture zone from the regional field is critical, in interpreting the 
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residual SEASAT geoid of Figure 3 chis is much loss cricical. This is so 
because the principal cause of the anomalies hounded by fracture ?.ones is the 
regional dlCturence in plate age. Each section ol plate liecu'een li- tj major t racciiLe 
zones is large enough ( = 1000 km) to be constantly near isostatic equilibrium. And 
each section of plate is large enough such that it has a characteristic 
geoid anomaly that, for the most part, depends critically on the absolute 
age of the plate. These points have been treated in some detail by Haxby 
and Turcotte (1978), Turcotte and McAdoo (1979), and Sandwell and Schubert 
(1980). Under the assumption of isostatic equilibrium, the geoid or gravity 
anomaly can be calculated from knowledge of the local vertical distribution 
of density as long as the depth of compensation is much less than the wave- 
length of the nomaly itself. For the plate regions between major fracture 
zones, the anomaly wavelength is ~1000 km. The geoid anomaly N(x,y) due 
to a density distribution p(x,y,z) over a depth L is given by (e.g. Turcotte 
and McAdoo, 1979) 

L • 

N = - J zAp(x,y,z)dz (!' 

o 

where z is the vertical coordinate, g is gravity at the surface and G is 
the universal gravitational constant. The application of (1) here is only an 

approximation'. To. model the details of the anomaly near fracture zones them- 
selves the more exact techniques of, say, Sandwell and Schubert (1982) can be 
employed. For the broad scale anomalies considered here, however, results 
gained from (1) and also empirical relations (see below) should suffice. 

i 

Once a function for p is known. throughout L, it is routine to compute 
N. For the present consideration, it is convenient to choose L as the 
thickness of the lithosphere. Since the thermal regime of the lithosphere 
changes with time, any equation of state will of necessity also be a 
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function of time. The result is not especially sensltve to the equation 
of state or to the assumed thermal dlstr Ihutlon for timen less than about 
50 m,y. 

For a thermal model of the lithosphere based on a thermal boundary layer, 
the geoid height relative to the ocean ridges is (Turcotte and McAdoo, 1979> 


N *= 


l.l7Gp„a(T,-T )L^ 

/ L O 

g 


1 + 


I-(PI_-P„> J 


( 2 ) 


where p ,is the mantle density at a depth L at the base of the lithosphere 

li 

where the temperature is is the surface temperature, a is the 

isothermal coefficent of thermal expansion, p^^ is the density of sea water 
and g is gravitational acceleration of the surface. A corresponding 
expression for the plate model is given by Sandwell and Schubert (1980) . 

2 

It is clear from (2) that the geoid anomaly is directly proportional to L , 

the thickness of the lithosphere, which for the thermal boundary 
2 2 

layer model is L = (2,32) Kt, where K is thermal diffusivity and t is time. 
With this substitution, (2) shows that the geoid anomaly is linearly dependent 
on Che absolute age (t) of the lithosphere. Here we are interested in the 
anomaly caused because of the Relative difference in age (At) of adjacent 
regions of plate bounded by fracture zones with large offsets (i.e, tens or 
more of millions o'f years) . This relative geoid anomaly (AN) is thus given 
by 


AN = - 6.30 


GBK 

g 


1 + 


2B 

TT(PL-f\,) 


At 


(3) 


where B * PrCifT ~T ). 

L L o 


The last term in brackets is always of order one and 
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U IS clenr l^hat AM OHK At/j;, which for 


) . 30 , <’ » 980 cm/s^ 


T -T = 1250°C,a» 3.3 x 10“^ flag’’^ and K. « 10 
L o 


cm'^/s.jAN j == 1.75 m per 


10 m.y. of offset. This time difference (At - 10 m.y.) is reasonable for 
many fracture zones and the anomaly is of the same order as those shovm by 


Figure 3. 

Although this result may apply when the lithosphere is young, 
an unattractive feature of it is the fact that the relative geoid anomaly 
is independent of the absolute age of the lithosphere itself, l^^ereas it 
has been shown by Sandwell and Schubert (1980) that the rate of change of 
the geoid anomaly across fracture zones decreases systematically with plate 
age. Of the analytical models, they found that the plate model matches the 
data better than the boundary layer model. But overall a good estimate of 
the geoid anomaly can be obtain'-d from the observed anomaly decay. Their 
actual data on the decay of the slope of the geoid anomaly with age for the 
North Atlantic is of the form (see their Figure 4) 


” « 2 X lO'"^ t - 0.2 
dt 


(4) 


where N is the geoid height in meters and t is age in millions of years. 
Integrating (4) and employing the condition that N(t=0) = 0 and subtx*acting 
two such formulas for adjacent regions of, respectively, age tj^ and 
gives a result for the relative geoid .anomaly. 


AN = (tj - t‘ 2 ) X 10 "^ - 0.2(t^ - t^) 


(5) 
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Beginning at a ridge with a relative offset of lOm.y. gives an anomaly 
of 1.9 m, which decays to about I m when the older place has an age of 60m.y.,and 
Che anomaly vanishes at an age of about 100 m.y. (The actual sign of the anomaly 
depends on the relative plate age.) F?om this result it is clear that geold 
anomalies related to regional changes in plate age will only he significant over 
places having an absolute age of less Chan about 80 m.y. The largest anomalies will 
of course be closest Co the ridge and they v.'ill have decayed to near nothing 
over SOro.y. In the Pacific significant anomalies can only be expected approxi- 
mately east of the longitude of Hawaii. Overall it is useful, however, to apply 
these results to the interpretation of the aerially large positive anomalies 
that lie near Che east Pacific rise; namely, Che anomaly northeast of Hawaii, 
the eastern end of the Line Islands anomaly, and the Pitcairn anomaly. 

The positive anomaly northeast of Hawaii is bounded by the Murray and Molokai 
fracture zones. Across the fracture zones in this vicinity the average offset in 
age is about 18 m.y. (Pitman et at., 1974). This would produce an anomaly of 
about 3 ra when the actual anomaly approaches 6 m. This is not necessarily 
a serious dissagreement, because the anomaly directly to the north, over the 
older Murray-Mendocino area, vt also positive and the relative difference is 
about 3 m. Why this positive spills northward to the Mendocino fracture zone 
and beyond may be due to Che proximity of North America. Nevertheless it seems 
that this relative. 3 m anomaly decreases markedly toward Hawaii. This may be 
due in part to a change in Che offset age as has been noted for Che Mendocino 
fracture zone by Sandwell and Schubert (1982). 

Moving south, the next section of plate with a major offset is bounded 
by the Molokai and Clipperton fracture zones. This place is relatively 
older and the geoid anomaly is slightly negative. The next positive anomaly 
is immediately to the south and is bounded, although not tightly, by the 
Clipperton and Galapagos fracture zones. The anomaly itself seems too broad 
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to be wholly the result of this relatively siimll offset ('^5 m.y.), wlilch 
would produce about 1 m anomaly when the actual anomaly has an average 
amplitude of about 2 m. The auuiMulv nuarly reaches the Mar«iuesas fracture, 
the next major fracture zone to the south. Moving southerly, normal to 
the fracture zones along about 235*^, each offset brings successively 
younger plate further west. The eastern portion of this Line Islands anomaly is 
probably the result of this cummulative effect it' offset. This anomaly also 
diminishes westward and then increases again, suggesting a change in the age of off- 
set, which along the Galapagos fracture zone may even change sign (Pitman et al, , 1974) . 

Pitcairn anomaly is the only one of these positive anomalies that is bounded 
by the east Pacific rise itself. It is bounded by the Marquesas fracture 
zone on the north and the Challenger on the south. The offset in age across 
these fracture zones is largest nearer the ridge and decreases westward. This 
has the effect of producing an anomaly nearer the ridge that attenuates 
rapidly to the west, The anomaly itself lies almost wholly on sea floor formed 
upon reorganization of the east Pacific ‘rise beginning' about 60m.y.b.p. (Herron, 
1972) . At this time the northwesterly trending ridge migrated north of the Eltanin 
fracture zone forming a new, northeast trending ridge and reorienting fracture 
zones to strike west-northwest. These newer fracture zones and young lithosphere 
produce the Pitcairn anomaly. 

Although these, interpretations in light of the presence of fracture 
zones having significant offsets in age seem reasonable, they are only 
qualitative. But they can be put on a firmer basis by computing a synthetic 
geoid based strictly on the relative ages of the Pacific plate; which is now in 
progress (Marsh and Hinojosa, 1983). .This is quite feasible for the Pacific 
plate because of the relative paucity of voluminous off-rldge volcanism 
for large distances near the west flank of the east Pacific rise. 
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Possible Size and Ascent Velocity of Hot Spot P, Indies 

Although tlie sncfcLcial expri'S.sion of !u»t spats ns voLf.nUsm and their 

kinematic implications are by now clearly evident, there is much less 
certainty about their subsurface structure and ultimate origin. That is 

little is known of the fluid mechanical conditions within the mantle that 
gives rise to these features. Morgan (1971, 1973) has suggested that they 
originate deep within the lower mantle as buoyant plumes, perhaps emanating 
from the region of the core-mantle boundary itself. Judging from the 
wide nature of lithospheric swells (»1000 km), and using a "simple aspect- 
ratio arguement", Crough (1978) has also suggested that these plumes originate 
from the lov;er mantle or at the core-mantle boundary itself. Sandwell 
(1982) has attempted to place bounds on the diameter of the hot-spot below 
the lithosphere by computing the gravity, geo id, topography, heat flow, 
and subsurface temperature as a function of diameter and plate velocity 
using an imposed temperature distribution. These results are not particu- 
larly diagnostic, but they suggest a sub-lithospheric hot-spot diameter of 
400-4000 km. It is the intent of this section to show some additional 
relationships between plume size, mantle viscosity, and heat transfer during 
ascent. 

A necessary assumption in relating hot-spots to a source region for 
plumes is that the plume must travel fast enough to the earth’s upper most 
regions without cooling and loosing its buoyancy. The rate of cooling 
depends critically on the diameter of the plume which is controlled by the 
viscosity field in the source region. Plumes rise as a result of a density 
inversion or gravitational instability. It has been shown through analysis 
and experiment that the characteristic size (a = radius) of a plume (of 
viscosity ^2) emanating from such a low density layer (of viscosity p^^) is 
given by (Marsh , 1979) 
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where is Che inicial thickness of the source layer. If « ^ 2 * the 
plume has a diameter about equal to the thickness of its source layer, 
which is approximately observed for salt domes within the earth's crust. 

A buoyant layer within the mantle could be produced by anomalous 
heating. The viscosity, which is othen^^ise nearly uniform (Cathles, 1973), 
could also be lowered by this heating. The reduction in viscosity, hov/ever, 
is unlikely to be by more than a factor of 100 (poise) (e.g. Harsh, 1982) 
unless it is partially melted, but there is no evidence that the 
temperature of the deep mantle is near its solidus. From Figure 5, 
computed from (4), for plume diameter is about three times 

the thickness of its source layer. Since a source layer much thicker than 

to 

about 100 km would be readily detectable by seismic methods, the plume 
diameter would be about 300 km; even for a viscosity contrast of 10^, the 
diameter is only about 400 km. The main lesson of equation (4) is that the 
plume diameter will be of the same order as the thickness of its source. 

Unless the source layer thickness is greater than about 5-10 times its 
depth, there are no results associated with the fluid mechanics of gravita- 
tional instabilities that give any information on the actual depth of the 
source (Marsh, 1975). The only region within the mantle where it is 
dynamically reasonable to suspect plume growth is within hot thermal 
boundary layers. Such a layer may occur at the core-mantle boundary 
(Elsasser et al., 1978), '^hich may be associated with the D" seismic layer. 
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although this is controversial. If the mantle has two or more tiers of 
convection cells, hot thermal boundary layers could also occur within the 
body of the mantle itself, but their existence is highly controversial. 

A small plume ascending from great depth may cool and loose its 
buoyancy before reaching the base of the lithosphere. Its final temperature 

depends critically on its rate of ascent, and even a small plume ascending suffi- 

« 

ciently fast will hardly cool at all. This relationship between ascent velocity, 

final temperature, and plume velocity may be investigated through a model of 

heat transfer. These models and methods have been extensively developed by 

Harsh (1978; 1982) and Marsh and Kantha (1978), for understanding 

the heat transfer of ascending bodies of magma. The method is parametric 

and general and may be used for any geometry, and only a brief description 

is given here. 

By conservation of energy, the mean temperature (T) In a plume changes 
with time in response to the total 'flux of heat (Q^) from the body and that due to 
adiabatic decompression (-yX). 

pCpV’ ^ « -Q^ r-Y(T)/ pCpV' (5) 


where p is density, Cp is specific heat, V' is volume and y is the usual adiabatic 
coefficient; all of which are considered constant (although this is not a necessary 
assumption). A dimensionless number involving can be conveniently defined as 



( 6 ) 


where Nu is the well known Nusselt number, a pure number, *and is the 
heat flux due to conduction (i.e. when the body is motionless). The flux 
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Q , AKcCT-T (t))/Ii, where A ir> the r.iirface area Involved in the lieat frouHt'er, 
t‘ci m 

Kc is the.rnial concluccivity , is the characteristic mantle geotherm iar 

from the plume, and L is the characteristic length scale of heat con- 
duction; for a plume L » radius (a). Substituting and (6) into (5) 
and rearranging gives (Marsh, 19 ?8) 

+ JT JT (t) -yT (7) 

at ra 

where J « Nu(AKc) /(LpCpV ) . This has the general solution 

T = I (8) 

J 

where J' = J + y» when coupled with an intial condition this describes the 
mean temperature of the body as a function of time or distance from its 
source. Since the final temperature must be greater then T^(L), where L is 
the total ascent distance, this places a condition on the ascent velocity. 
Although the ascent velocity has not yet appeared explicitly, it enters 
through the relationship for Jt. 

For magmatic transfer through the lithosphere, solutions to (8) are 
given by Marsh (1978; 1982) and Marsh and Kantha (1978) where transfer by 
diapirism, stoping, zone melting, and pipe flow are considered. Here we 
consider an anomalous hot plume ascending from a buoyant region. The region 
itself need not be globally continuous , but only locally con- 
tinuous. Its mean temperature is Tp, which is the temperature of the plume 
itself at its source. The normal mantle temperature at this depth is T^. 
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The mantle temperatuire ts taken lo he adlnhatle throuftlunit Che mantle, and 
hence varies as l'^(t)/T^ -- with depth (or time) until near the 

uppermost niontie. Substitucinij into (8) gives upon incegratton 

^ « EXP(-(Yt^)t')-EXP(~((J-H)t^)tO'<-(T /t^)EXF(-((J+Y)t^)t’). (9) 

o ^ 

where the intial condition T(0) « Tp has been used. The mean temperature 

of the plume as a functionof t' ( » t/t^) and Jt^ as described by (9) is 

given by Figure 6. Both the mantle and the plume have been assumed to 

have the same adiabatic coefficient (y = agt^/Cp ® 0.7, where Cp is 

specific heat and other symbols as before) . 

The results of Figure 6 suggest that if the plume is to remain hot 

and buoyant and reach the upper mantJn Jt^ = 1, from which the ascent 

velocity can be calculated. Recalling that J 3 (NuAKc)/(LpCpV’) , where 

2 

for a plume A/V = 2TraA/ira Z = 2/a , where a is radius and i. is a length, 

-22 -22 
Kc/pCp - 10 cm /s, and L = a, we have J = Nu2xl0 /a . It is clear that 

Nu must now be used and several choices are possible. The Nu for a plume 

can be approximated by that for flow in a pipe, which for a fully developed 

pipe flow is well known to be 4.36 (e.g. Kays, 1966). TMs result holds 

only far from the thermal entry region of the pipe. The length of the 

thermal entry region for most geophysical flows is large (Marsh, 1978) and 

2 —7 7 

it is proportional to Va /K, where V is velocity. If V = 10 cm/s, a = 10 cm 

-2 2 4 

and K - 10 cm /s, the thermal entry length is about 10 km, and for larger 

velocities it increases. Thus Nu = '4.36 is a conservative approxiamtion for heat 

transfer, more chan likely it will be significantly larger, even by a factor of ten. 
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Then, .1 B.7?.*<10 /.i' mul fhe tut. si cinsn Lh j'iven l»y lI5H;s , 

where H Is the mimertcal value of Jt from tho curves oi‘ Fiuure 6. 

o 
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For Jt^«l«H and, from Figure 5 for, say, n“200 km, the ascent time fs t^«4xl0 sec, 
which for an ascent distaiu.' of 2500 km, gives a typical ascent velocity of 
about 2 cm/yr. That is, if a plume of a diameter of 400 kin ascends a distance 
of 2500 krti at a velocity of about 2 cm/yr, it will still be anomalously 
hoc when it arrives in the upper mantle. It must be reminded, however, that 
this is an absolute mimimum velocity. For a more reasonable Nu (~10x larger) 
the velocity must be 20 cm/yr, and if the radius is 100 km, V «80 cm/yr. 

The main result of this calculation is chat it is apparently possible even for 
relatively small plumes to ascend a large distance through the mantle without 
totally loosing their original anomalous temperature. It should also be 
pointed out, however, that even an adiabatic plume cools more than an equivalent 

"normal" mantle during its ascent. This is because the adiabatic temperature 
change is proportional to the temperature of the body itself. Hence even for 
identical thermal properties, a hotter body cools faster than an equivalent 
cooler body. 


I. Conclusions 

The gravity field over the Pacific determined by Sjogren's method with 
about 90 passes of CEOS-3/ATS~6 satellite-to-satelllCe tracking data. This 
new determination is in good agreement with our previous (1981) determination, 
of the gravity field in this area using 40 passes of SST data. A comparison 
of this map with the more conventially-determined GEM satellite gravity field 
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.‘Utnwa fjood nr,r(»etnpnK We nl{»o ohow Llie Meold over llu* I’iU-lfio th*f;ermln«tl 

Uoiii i<KASA'I alUiuet.er data and it too auteea well with both OHM and SST, 

but chia ficoid contains far more detail than either of these other maps. 

It has been previously noted that areas with residual depth anomalies 

associated with hot-spots correlate with gooid and gravity anomalies. These 

so-called swells do often correlate with the present maps, but there are 

some clear exceptions and complications. The Marquesas-line swell of Crough 

and Jarrard (1981) was found by them to correlate well with the CEOS-3 geoid 

of Brace (1977). The positive geoid anomaly near the Marquesas Islands 

of that work, however, was not found in our earlier SST study (Marsh et al., 

1981) and the present fields verify its absence. This errant positive 

anomaly is apparently the result of poor GEOS-3 coverage in this area. 

Those anomalies not associated with thermal swells correlate well with 

relatively young areas of the seafloor bounded by fracture zones. Because 

the major fracture zones of the Pacific are fairly evenly spaced and trend 

approximately east-west, they produce a similar fabric in the gravity and geoid 

fields. The anomalies are essentially framed by these fracture zones. For 

a typical offset in age of -15 m.y. the resulting geoid anomaly is about 2.5 m, 

which is close to that observed near the east Pacific rise. Because the traces 

of hot spots are not parallel to fracture zones, there is an interference in 

» 

anomalies from each source. The resulting anomaly field may thus be separable 
using the known ages and history of the Pacific plate. Early modeling shoves 
that because the anomalies due to age offsets attenuate with absolute age, 
large geoid anomalies over older sea floor ('v70 m.y. like NE of Hawaii) can not 
be completely explained by this offset alone (Marsh and Hinojosa, 1983). A part 
of the cause may also be below the lithosphere. 
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An oxai'ilnnl Ion of tlin nixr* .uul n'lftMif volo«-(tv noenKnnry to briny, mantli? 
plniiiwii to the upptu* iiiantli* wlthonr (Hujlinji fjliowo titat pluiiit* diai.wtwr in apt in 
be of the name order ns the chlcknens of Its .source. Because lou density 
reglon.s within the mantle thicker than about 1000 km, were they to exist, .should 
have been discovered through seismology, plume diameter is probably limited to less 
than about 400 km. And the ascent velocity needed to prevent complete cooling 
Is at least 3 cm/yr; it could be ten times larger. 

The explanation of this distinctive pattern of gravity and geoid 

anomalies in the Pacific has been of interest since the first indication of 

its presence (Marsh and Marsh, 1976). At that time we suggested they 

might reflect the presence of a small scale form of convection in the uppermost 

> 

mantle. Because several of these anomalies begin very near the east Pacific 
rise and because some time is necessary to initiate small scale convective rolls, 
McKenzie et al. (1980) suggested that these anomalies must arise from small 
scale convective instabilities in the lo^Aer limb of a larger cell confined to 
the upper mantle. These secondary instabilities would then be already estab- 
lished as they reached the ridge itself. The close correlation between 
sea-floor age, fracture zones and these near-ridge anomalies makes the inter- 
pretation of McKenzie et al. (1980) untenable. In fact as it now appears 
the structure of the lithosphere itself may be a major factor in causing 
these distinctive anomalies. Nevertheless, the ultimate origin of both hot 
spots, the regularity of fracture zones with major offsets, and some geoid 
anomalies over old, smooth sea floor is still unknown, and they could 
conceivably be related to the effect.? of small convection. 
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Fi{»ure Captions 

. Distribution of SST tracks used in this study. 

. SST gravity map for the central Pacific region. The circles 
have radii of 30“and 40° about the subsatellite point, and witliin 
this range the gravity is essential radial. 

. SEASAT altimeter geoid for the Pacific relative to the twelfth 
degree and order GEM lOB geoid. The trace of the east Pacific 
rise is also shown (far right) as well as the 20 and 65 m.y. 
isochrons for the age of the Pacific plate. 



The SEASAT geoid of Figure 3 superimposed on the map of the age 
of the ocean basins of Pitman et al., (1974). 

The relationship between the diameter (2a) of a plume rising from 
a low density region of thickness and viscosity p 2 relative to 
a surrounding imantle of viscosity The right axis gives diameter 

for a source thickness of 100 km. 

The mean temperature (T/T^). of a section of plume as it ascends 
through mantle whose temperature is adiabatic and described by 
Tm(t). Only for values of Jt^ = 1 does the plume reach the upper 


mantle still anomalously hot. 
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Figure 2. SST gravity map for the central Pacific region. The circles have 
radii of 30* and 40* about the subsatellite point, and within thi 
range the gravity is essential radial. 





SEASAT MEAN SEA SURFACE 
MINUS GEM 10B (12,12) GEOID 
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PLUME DIAMETER 



VICOSITY CONTRAST 


Figure 5. The relationship between the diameter (2a) of a plume 
rising from a lot* density region of thickness h2 and 
viscosity P2 relative to a sourroundlng mantle of 
viscosity uj. The right axis gives diameter for a 
source thickness of 100 km. 


diameter (hglOO km) 



PLUME TEMPERATURE 



Figure 6. The mean temperature (T/Tq) of a section of plume as It ascends through nant 
whose temperature Is adiabatic and described by Tm(t). Only for values of 
Jt =* I does the plume reach the upper mantle still anomalously hot. 




